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Thermodynamic  properties  and  detailed  chemical  kinetic  models  have  been  developed  for  the  combus¬ 
tion  of  two  oxygenates:  methyl  butanoate,  a  model  compound  for  biodiesel  fuels,  and  methyl  formate,  a 
related  simpler  molecule.  Rond  additivity  methods  and  rules  for  estimating  kinetic  parameters  were 
adopted  from  hydrocarbon  combustion  and  extended.  The  resulting  mechanisms  have  been  tested  against 
the  limited  combustion  data  available  in  the  literature,  which  was  obtained  at  low  temperature,  subat- 
mospheric  conditions  in  closed  vessels,  using  pressure  measurements  as  the  main  diagnostic.  Some  qual¬ 
itative  agreement  was  obtained,  but  the  experimental  data  consistently  indicated  lower  overall  reactivities 
than  the  model,  differing  by  factors  of  10  to  50.  This  discrepancy,  which  occurs  for  species  with  well- 
established  kinetic  mechanisms  as  well  as  for  methyl  esters,  is  tentatively  ascribed  to  the  presence  of  wall 
reactions  in  the  experiments.  The  model  predicts  a  region  of  weak  or  negative  dependence  of  overall 
reaction  rate  on  temperature  for  each  methyl  ester.  Examination  of  the  reaction  fluxes  provides  an  expla¬ 
nation  of  this  behavior,  involving  a  temperature-dependent  competition  between  chain-propagating  uni- 
molecular  decomposition  processes  and  chain-branching  processes,  similar  to  that  accepted  for  hydrocar¬ 
bons.  There  is  an  urgent  need  to  obtain  more  complete  experimental  data  under  well-characterized 
conditions  for  thorough  testing  of  the  model. 

Introduction  fuels  consist  of  mixtures  of  saturated  and  unsatu- 

,7  1  1  .  1  r  t  ui  -1  ]  ■  I  r  4.  rated  methyl  esters  containing  carbon  chains  12  or 

F  uels  derived  from  vegetable  oils  and  animal  tats  ^  i  m 

offer  an  alternative  to  petroleum  products  for  diesel  [1]'  Numerous  engine  corn- 

engine  combustion.  These  biodiesel  fuels  are  attrac-  studies  have  been  performed  with  these  fu- 

tive  because  of  their  low  sulfur  content,  which  allows  *eir  constituents  (e.g.,  [3]).  However,  the 

the  use  of  a  catalyst  to  remove  NO,  from  lean-burn  few  fundamental  studies  suitable  for  comparison 
engine  exhaust  to  meet  new  stringent,  governmental  chemical  kinetic  models  are  restricted  to  the 

emission  standards.  Biodiesel  fuels  are  also  renew-  pyrolysis  regime  (e.g.,  [4]),  as  are  the  only  known 
able  fuels  which  have  low  impact  on  global  warming  detailed  chemical  kinetic  models  (e.g.,  [4]). 
and  can  help  limit  dependence  on  foreign-derived  This  paper  describes  the  development  of  a  de¬ 
fuel  supplies.  With  oxygen  content  typically  10%  or  tailed  chemical  kinetic  model  for  methyl  butanoate 

greater  by  mass  [1],  biodiesel  fuels  may  also  provide  (nC3H7C(  =  0)0CH3),  a  surrogate  for  biodiesel  fu- 

soot-reduction  benefits  similar  to  those  observed  for  els.  Although  methyl  butanoate  does  not  have  the 

other  oxygenated  diesel  fuels  and  additives  [2].  Al-  high  molecular  weight  of  a  biodiesel  fuel,  it  has  the 

though  oxygenated  fuels  do  not  directly  reduce  NO,  essential  chemical  structural  features,  namely  the 

emissions,  their  soot  reduction  capabilities  increase  RC(  =  OlOCHg  structure  (where  R  is  an  alkyl  or  al- 

flexibility  in  the  choice  of  engine  operating  condi-  kenyl  radical).  Little  literature  information  on  ther- 

tions  for  controlling  NO,  emissions.  modynamics  and  kinetics  of  RC(  =  0)0CH3  and 

For  engine  applications,  biologically  derived  fatty  species  derived  from  it  is  available.  This  paper  strives 
acids  are  ordinarily  converted  into  methyl  or  ethyl  gQ  thermodynamic  and  chemical 

esters  to  improve  physical  properties  such  as  viscos-  gi^ggg  Methyl  butanoate  was  chosen  as  a  sur- 
ity  and  melting  and  boiling  points.  Typical  biodiesel  ^^g^jg  ^^gigg^ig  Hrger  methyl  esters  in  order  to 

_  obtain  a  reaction  mechanism  of  manageable  size. 

*Permanent  address:  Chemistry  Department,  Galway-  Yet,  methyl  butanoate  is  large  enough  to  allow  the 
Mayo  Institute  of  Technology,  Galway,  Ireland.  fast  RO2  isomerization  reactions  important  in  the 
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TABLE  1 

Hydrogen  Bond  Energies  and  Increments 


Group  Name 

^Hbe(298)“ 

(kcal/mol) 

^5(298)*’ 
(cal/(mol  K)) 

Basis  for  Values 

CH30C(-)  =  0 

98.4 

0.84 

Comparison  of  HOCO  and  HOCHO  BAC-MP4 
calculations  [20] 

•CH2OCHO 

100.35 

1.73 

AH:  CBS-Q  calculations  with  group  balance 
isodesmic  reactions  [18];  AS:  •CH2OR  group  [16] 

•CH2C(  =  0)0R 

98.8 

-1.31 

Comparison  of  CH2COOH  and  CH3COOH  BAC-MP4 
calculations  [20] 

CH3CH(-)C(  =  0)0R 

96.2 

-2.37 

Values  for  -CHaC]  =  0)OR,  changed  by  the  difference 
between  CH3CH(-)CHO  and  -CHaCHO  [16] 

“Bond  energy,  for  example,  CHsCHl  =  O)  =>  CHsOC-  =  O  +  H,  dff 2^(298)  =  98.4. 

'’Increment:  value  added  to  parent  (Rln2  for  electron  spin  and  symmetry  corrections  are  made  on  a  molecule-specific 
case  in  the  THERM  program  [16]). 


low-temperature  chemistry  that  controls  fuel  autoig¬ 
nition  under  conditions  found  in  diesel  engines.  A 
model  for  methyl  formate  (HC(  =  0)0CH3)  devel¬ 
oped  using  the  same  reaction  rate  rules  is  presented 
as  well. 

Both  methyl  butanoate  and  methyl  formate  are 
the  subject  of  several  low-temperature  oxidation 
studies  performed  in  small,  constant- volume,  iso¬ 
thermal  static  reactors,  with  O2  as  the  oxidant  [5-8], 
Temperatures  and  initial  pressures  range  from  520 
to  740  K  and  from  13  to  54  kPa,  respectively.  Most 
experiments  use  rich  mixtures,  but  some  experi¬ 
ments  are  performed  with  the  equivalence  ratio  as 
low  as  0.65  (methyl  butanoate)  and  0.1  (methyl  for¬ 
mate).  Pressure  is  usually  the  only  experimental  di¬ 
agnostic,  with  results  presented  in  the  form  of  in¬ 
duction  time  or  maximum  rate  of  pressure  rise, 
taken  as  representing  the  maximum  reaction  rate. 
Limited  major  species  measurements  [5,7,9]  are 
available  for  methyl  formate. 

Methyl  butanoate  [6,8]  exhibited  a  negative  tem¬ 
perature  coefficient  (NTC)  region,  with  a  peak  in 
maximum  reaction  rate  occurring  at  about  585  K. 
No  NTC  behavior  was  observed  for  shorter-chain 
methyl  esters  including  methyl  formate  [6].  Associ¬ 
ated  with  NTC  behavior  is  the  occurrence  of  cool 
flames  [10].  Cool  flame  limits,  based  on  observations 
of  luminosity,  were  reported  for  methyl  butanoate 
[8,11],  while  no  cool  flames  were  seen  for  shorter- 
chain  methyl  esters  [8]  including  methyl  formate 
[11].  More  recently  Baronnet  and  Brocard  [12] 
summarized  cool  flame  and  NTC  observations  and 
proposed  a  criterion  for  their  occurrence  based  on 
the  availability  of  relatively  fast  isomerization  path¬ 
ways  for  alkyl  peroxy  radicals  derived  from  the  fuel. 
Since  several  alkyl  peroxy  radicals  derived  from 
methyl  butanoate  can  participate  in  rapid  isomeri¬ 
zation  reactions,  the  criterion  successfully  predicts 
NTC  behavior  and  cool  flames  for  this  compound. 


However,  methyl  formate  meets  the  criterion  as  well 
through  transfer  of  an  aldehydic  hydrogen  atom  via 
a  6-member  ring.  Baronnet  and  Brocard  attribute  its 
lack  of  NTC  and  cool  flame  behavior  to  its  low  en¬ 
thalpy  of  combustion.  In  recent  work.  Good  and 
Francisco  [13]  investigated  the  oxidation  mechanism 
of  methyl  formate  using  ab  initio  computational 
chemistry  methods. 


Model  Development 

Chemical  kinetic  calculations  were  performed 
with  the  SENKIN  code  [14],  which  requires  ther¬ 
modynamic  and  kinetic  data  as  input.  Thermody¬ 
namic  properties  of  the  methyl  esters  and  their  de¬ 
composition  products  were  estimated  using  group 
additivity  [15]  as  implemented  in  the  THERM  code 
[16]  with  updated  group  values  [17-19].  Several  new 
bond  dissociation  groups  developed  to  predict  radi¬ 
cal  properties  are  listed  in  Table  1.  During  the  de¬ 
velopment  of  the  methyl  butanoate  mechanism,  the 
value  of  the  C-H  bond  dissociation  energy  for  the  C 
adjacent  to  the  carbonyl  was  originally  taken  from 
Bozzelli  and  Chen  [18]  (93.1  kcal/mol).  However, 
the  resultant  mechanism  predicted  no  NTC  behav¬ 
ior  for  methyl  butanoate.  In  the  final  version  of  the 
mechanism,  the  C-H  bond  strength  was  increased 
to  96.2  kcal/mol,  a  value  derived  using  BAC  MP4 
estimates  [20]  of  similar  radical  species  (Table  1). 
This  increase  in  bond  strength,  with  all  the  required 
accompanying  changes  in  reaction  rate  constants,  re¬ 
produced  the  experimentally  observed  NTC  behav¬ 
ior. 

The  full  chemical  kinetic  mechanism  is  available 
from  E.  M.  Fisher  (email:  emf4@cornell.edu).  For 
the  most  part,  it  follows  the  rules  developed  by  Cur¬ 
ran  et  al.  [21]  for  n-heptane,  as  updated  in  an  iso¬ 
octane  study  [19].  Because  the  C-H  bond  strength 
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Fig.  1.  Calculated  values  (lines)  and  experimental  [6] 
values  (symbols)  of  dP/dt^^  for  isothermal,  constant-vol¬ 
ume  combustion:  1:2  mixture  of  methyl  butanoate  and  O2; 
initial  pressure  40.0  kPa.  1:1  mixture  of  methyl  formate  and 
O2;  initial  pressure  53.3  kPa. 


for  the  C  adjacent  to  the  carbonyl  in  methyl  butan¬ 
oate  was  very  close  to  the  value  for  tertiary  carbon 
(96.2  vs.  96.5  kcal/ mol),  the  rules  for  tertiary  carbon 
were  used  for  this  site.  The  following  classes  of  re¬ 
action  were  added  or  changed:  (1)  Fuel  H  abstrac¬ 
tion  by  RO2  species  formed  from  the  fuel  was  given 
half  the  rate  for  H  abstraction  by  HO2.  (2)  Abstrac¬ 
tion  of  the  aldehydic  H  in  species  with  two  carbonyl 
groups  was  given  rates  based  on  aldehydic  H  ab¬ 
straction.  (3)  A  new  approach  for  OH  and  O  addition 
to  olefinic  species  was  based  on  propene  literature 
[22-24].  (4)  For  radical  attack  on  cyclic  compounds, 
Evans-Polyani  relations  were  used  for  activation  en¬ 
ergy  estimates,  representing  each  C-H  bond  94  kcal/ 
mol  or  weaker  with  its  own  set  of  abstraction  reac¬ 
tions.  In  some  cases,  the  product  distribution  was 
chosen  for  convenience  rather  than  to  represent 
products  specific  to  the  radical  attack  site. 

The  following  changes  were  made  in  the  rules  for 
obtaining  rate  constants:  (1)  For  isomerization  of  es¬ 
ter  radicals,  values  of  ring  strain  contributions  to  ac¬ 
tivation  energy  recommended  by  Benson  for  alkyl 
radicals  [15]  were  used;  for  the  RO2  (alkyl  peroxy 
analog)  and  O2QOOH  (peroxy-hydroperoxy-alkyl 
analog)  species  derived  from  the  esters,  Currans 
ring  strain  values  were  used  [21].  (2)  Preexponential 
factors  for  isomerization  of  these  three  classes  of 
species  were  obtained  respectively  by  RADICALC 
[25]  calculations  of  the  transition  state  by  matching 
rate  constants  to  those  of  Curran  et  al.  at  700  K  and 
by  adjusting  preexponential  factors  for  RO2  isom¬ 
erization  to  account  for  new  degeneracy.  (3)  For 
QOOH  (hydroperoxyalkyl  analog)  scission  to  form 
olefin,  carbonyl,  and  OH,  a  reverse  reaction  rate 
comparable  to  that  for  the  reverse  of  alkyl  radical 
beta  scission  was  used.  (4)  Several  specific  reaction 


rates  were  chosen  by  analogy  with  rates  for  oxygen¬ 
ated  species,  instead  of  using  the  rules  developed  for 
alkanes  and  their  derivatives;  these  reactions  fall 
mainly  into  the  following  categories:  beta  scission, 
unimolecular  fuel  decomposition,  and  RO  (alkoxy 
analog)  decomposition. 

The  methyl  butanoate  and  methyl  formate  sub¬ 
mechanisms  were  added  to  an  updated  version  [19] 
of  an  existing  n-heptane  mechanism  [21].  Reactions 
involving  larger  hydrocarbons  (C4  and  higher)  were 
removed  to  speed  calculations.  The  exclusion  of 
higher  hydrocarbons  had  a  minimal  effect  on  cal¬ 
culated  results  with  these  ester  fuels. 


Computational  Results:  Comparison  with 
Literature  Data 

Calculations  were  performed  for  comparison  with 
static  reactor  experimental  results  [5-8]  following 
the  idealized  description  of  experimental  conditions 
provided,  namely,  combustion  was  assumed  to  take 
place  isothermally  and  at  constant  volume.  Two 
types  of  results  were  compared:  the  maximum  value 
of  the  derivative  of  pressure  with  respect  to  time 
(dP/dt^^),  and  the  induction  time,  defined  as  the 
time  at  which  dP/dt^^  occurred  [6].  dP/dt^^,  which 
Parsons  and  Danby  [5]  indicated  is  related  to  the 
maximum  reaction  rate,  was  obtained  from  calcu¬ 
lated  pressure  by  numerical  differentiation.  In  the 
experiments,  pressure  was  measured  with  a  manom¬ 
eter  (probably  recorded  manually)  and  thus  the  data 
would  not  have  allowed  resolution  of  pressure 
changes  more  quickly  than  every  several  seconds; 
computed  pressure  derivatives  did  not  take  into  ac¬ 
count  the  manometer  response  time. 

Calculations  consistently  indicated  that  combus¬ 
tion  was  more  rapid  than  was  observed  experimen¬ 
tally.  Fig.  1  shows  the  calculated  and  experimental 
[6]  values  of  dP/dt^.^  for  the  two  methyl  esters,  as 
functions  of  the  inverse  of  temperature.  For  com¬ 
parison,  Fig.  2  shows  the  same  type  of  experimental 
data  for  propane  and  propene  [26],  along  with  com¬ 
putational  results  obtained  with  the  hydrocarbon 
submechanism  used  in  the  ester  calculations.  The 
present  reaction  mechanism  is  based  on  ones  used 
successfully  to  model  propene  oxidation  in  closed 
vessels  [24]  and  low-temperature  chemistry  of  pro¬ 
pane  in  flow  reactors  [27,28].  The  experimental  ap¬ 
paratus  and  methods  used  in  the  various  experi¬ 
ments  shown  in  Figs.  1  and  2,  performed  in  the  same 
research  group,  were  very  similar.  F or  all  four  com¬ 
pounds,  the  calculated  dP/dt^,^  was  higher  than  the 
experimental  value,  exceeding  it  by  factors  of 
roughly  10,  50,  10,  and  50,  respectively,  for  methyl 
butanoate,  methyl  formate,  propane,  and  propene. 
Features  of  the  experimental  curves  appear  to  be 
displaced  50  to  70  K  higher  in  the  calculations,  pos¬ 
sibly  indicating  that  heat  release  in  the  experiments 
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Fig.  2.  Calculated  values  (lines)  and  experimental  [26] 
values  (symbols)  of  dP/dt^^  for  isothermal,  constant-vol¬ 
ume  combustion.  3:25  mixture  of  propane  or  propene  with 
O2;  initial  pressure  37.3  kPa. 


Fig.  3.  Calculated  values  (lines)  and  experimental  [6] 
values  (symbols)  of  the  induetion  time  for  methyl  butan- 
oate  (a)  and  methyl  formate  (b).  Same  conditions  as  for 
Fig.  1. 


caused  deviations  from  the  isothermal  assumption 
made  in  the  model.  A  comparison  to  a  different 
static  reactor  experimental  data  set  [8]  gave  quali¬ 
tatively  similar  results:  calculated  values  of  dP/dt^^ 
showed  a  similar  shape  to  experimental  data,  but 
with  values  100  times  higher. 

Discrepancies  in  the  overall  speed  of  reaction,  as 
represented  by  dP/dt^^,  may  be  due  to  errors  in  the 
kinetic  mechanism  or  to  departures  of  the  experi¬ 
ment  from  ideal  conditions.  Because  propane  and 
propene,  which  have  relatively  well-established  ki¬ 
netic  mechanisms,  have  discrepancies  comparable  to 
those  for  the  methyl  esters,  a  problem  with  the  ex¬ 
perimental  methodology  is  likely.  One  plausible  ex¬ 
planation  is  the  presence  of  wall  reactions  in  the  ex¬ 
periments,  consistent  with  the  reported  increase  in 
reactivity  with  decreasing  surface-to-volume  ratio  in 
other  experiments  of  the  Mulcahy/Parsons  research 
group  [29].  Cheaney  et  al.  [30]  report  that  HO2  re¬ 
acts  on  the  surface  of  treated  silica  vessels  which  as 
a  result  of  treatment  have  thin  films  of  silicic,  boric, 
or  other  acids.  They  postulate  that  HO2  donates  a 
proton  to  the  surface.  It  is  likely  that  other  peroxy 
radicals  would  behave  in  a  manner  similar  to  HO2. 
Although  vessel  surface  preparation  was  not  dis¬ 
cussed  by  Parsons  et  al.  [5-7]  or  Mulcahy  [26], 
closed  vessels  are  often  treated  with  acids  prior  to 
their  experimental  use  (e.g.,  [24]).  The  magnitude 
of  the  effect  of  peroxy  radical  losses  was  investigated 
computationally  by  performing  additional  calcula¬ 
tions  with  a  rough  representation  of  this  type  of  wall 
reaction.  Reactions  converting  all  peroxy  and  per- 
oxy-hydroperoxyalkyl  radicals  into  fictitious  non-re¬ 
active  species  were  added  to  the  mechanism.  Rates 
for  these  reactions  were  chosen  to  be  equal  to  the 
collision  rate  of  the  species  on  the  wall,  multiplied 
by  a  constant  “efficiency.”  The  vessel  volume  (0.25 
L)  was  assumed  to  be  spherical  and  well  mixed.  An 
efficiency  of  1.3%  for  all  relevant  species  brought 
induction  time  and  dP/dt^^  into  agreement  with  ex¬ 
periment  near  the  middle  of  the  experimental  tem¬ 
perature  range  for  methyl  butanoate.  However,  it 
overcorrected  the  methyl  formate  dP/dt^^  results. 

Further  indications  of  overall  reactivity  are  pro¬ 
vided  in  Fig.  3,  which  shows  calculated  and  observed 
[6]  induction  times  for  the  two  methyl  esters.  Com¬ 
parable  experimental  data  is  not  available  for  pro¬ 
pane  and  propene.  Two  different  measures  of  in¬ 
duction  time  were  applied  to  the  calculated  data:  (1) 
the  time  at  which  the  maximum  value  of  dP/dt  oc¬ 
curs,  and  (2)  the  time  at  which  half  the  fuel  is  con¬ 
sumed.  These  two  measures  give  similar  answers  at 
high  and  low  temperatures  but  differ  at  intermediate 
values.  The  induction  time  based  on  the  fuel  de¬ 
struction  exhibits  a  smooth  dependence  on  tem¬ 
perature,  showing  an  NTC  or  near-NTC  region  in 
which  reactivity  has  a  low  or  negative  dependence 
on  temperature.  In  contrast,  the  data  based  on  dP/ 
dt  shows  an  unusual  dependence  on  temperature. 
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Fig.  4.  Calculated  values  (lines)  and  normalized  exper¬ 
imental  [6]  values  (symbols)  of  dP/dt^^  for  isothermal, 
constant-volume  combustion  at  O2;  pressure  26.7  kPa, 
temperature  653  K. 


Instead  of  a  region  of  NTC  behavior  as  temperature 
increases,  there  is  an  instantaneous  shift  to  a  longer 
induction  time.  The  shift  is  due  to  the  shape  of  the 
pressure  histories.  Both  methyl  esters  show  complex 
dP/dt  curves  as  functions  of  time,  with  multiple 
peaks  at  certain  temperatures.  The  relative  promi¬ 
nence  of  these  peaks  changes  as  temperature  in¬ 
creases,  leading  to  an  abrupt  change  in  induction 
time  as  the  location  of  the  maximum  value  of  dP/dt 
shifts  from  an  earlier  peak  to  a  later  one.  Because 
the  half  width  of  the  first  peak  is  so  short  (less  than 
a  second  typically),  it  is  likely  that  it  would  not  have 
been  observed  experimentally  with  a  manometer. 
The  other  measure  of  induction  time  appears  pref¬ 
erable  for  this  reason,  and  also  because  it  produces 
NTC  results  that  agree  with  those  seen  in  Fig.  1. 
Consistent  with  the  dP/dt^^  results,  calculated  in¬ 
duction  times  obtained  with  either  measure  are  fac¬ 
tors  of  10  to  50  times  shorter  than  those  observed 
experimentally.  For  both  esters,  calculated  induction 
times  with  wall  reactions  agreed  with  experiment  at 
some  point  in  the  experimental  range.  However,  the 
model  predicted  a  steeper  temperature  dependence 
of  both  induction  time  and  dP/dt^^  than  observed 
in  the  experiments. 

Qualitative  agreement  was  obtained  with  experi¬ 
ments  [7]  on  the  effect  of  ester  pressure  and  O2 
pressure  on  reactivity  at  constant  temperature.  Fig. 
4  shows  dP/dt^^  calculations  along  with  experimen¬ 
tal  data  that  has  been  normalized  to  match  calcula¬ 
tions  at  one  point  on  each  curve.  The  shape  of  the 
curves  for  both  methyl  esters  is  well  predicted.  Simi¬ 
lar  agreement  was  obtained  for  the  effect  of  ester 
pressure  on  induction  time  based  on  fuel  consump¬ 
tion.  For  the  effect  of  O2  pressure  on  dP/dt^,^, 
agreement  in  the  curve  shapes  was  good  for  methyl 
butanoate  but  poor  for  methyl  formate. 


Results:  NTC  Behavior 

Regions  either  of  NTC  behavior  or  of  a  related 
weak  dependence  of  overall  reaction  rate  on  tem¬ 
perature  (near-NTC  behavior)  can  be  identified  by 
examining  either  induction  time  or  dP/dt^^  (indic¬ 
ative  of  maximum  reaction  rate).  Calculated  dP/ 
dtm-jx  results  indicate  that  propane,  propene,  methyl 
butanoate,  and  methyl  formate  have  regions  of  NTC 
or  near-NTC  behavior  located  at  665-715,  690-765, 
640-670,  and  690-750  K,  respectively.  The  experi¬ 
mental  data  show  similar  behavior  for  propane  and 
methyl  butanoate  at  50  to  70  K  lower  temperatures. 
Propene  experimental  data  are  ambiguous,  showing 
an  inflection  that  can  be  intei'preted  as  near-NTC 
behavior.  Taking  into  account  the  temperature  shift 
between  experiment  and  calculation,  methyl  for¬ 
mate  experiments  appear  to  have  been  performed  at 
slightly  too  high  a  temperature  to  observe  the  mild 
near-NTC  region  seen  in  the  calculations.  Induction 
times  based  on  fuel  destruction  show  NTC  or  near- 
NTC  regions  consistent  with  those  of  Fig.  1,  while 
experimental  induction  times  show  no  NTC  behav¬ 
ior  over  their  temperature  range. 

Reaction  pathway  analysis  provides  an  explanation 
for  the  methyl  butanoates  NTC  behavior  similar  to 
that  accepted  for  hydrocarbon  fuels,  (e.g.,  [21]). 
NTC  behavior  can  be  explained  by  examining  the 
net  creation  and  destruction  of  radicals  in  pathways 
available  to  the  fuel  molecule  and  its  derivatives.  In 
this  isothermal  system,  chain  branching  processes 
(which  increase  the  radical  pool)  have  a  particularly 
important  effect  on  the  overall  reaction  rate.  At  two 
points  in  the  combustion  process,  there  is  a  com¬ 
petition  between  (1)  unimolecular  decomposition 
pathways  leading  almost  exclusively  to  chain  propa¬ 
gation  reactions,  and  (2)  addition  to  O2,  which  leads 
much  more  frequently  to  chain  branching  reactions. 
Fig.  5  shows  this  competition  for  the  radical  formed 
by  abstraction  of  one  of  the  H  atoms  beta  to  the 
carbonyl  group  in  methyl  butanoate,  one  of  the  four 
possible  radicals  of  methyl  butanoate.  The  first  com¬ 
petition  occurs  for  radicals  formed  directly  by  H  ab¬ 
straction  from  the  fuel  (reaction  1,  Fig.  5).  These 
species  can  either  decompose  by  beta  scission  (re¬ 
action  2)  to  produce  one  radical  and  one  stable  mol¬ 
ecule,  or  can  add  to  O2,  which  in  some  cases  (see 
below)  produces  more  than  one  radical  (reaction  3). 
The  second  competition  occurs  for  the  QOOH  spe¬ 
cies  formed  by  isomerization  (reaction  4)  of  the  O2 
addition  products  of  the  initial  radicals.  These  spe¬ 
cies,  too,  can  either  decompose  unimolecularly  (re¬ 
action  5)  to  a  radical  and  either  a  cyclic  ether  species 
or  an  olefinic  species,  or  can  add  to  O2  (reaction  6). 
The  products  of  this  second  O2  addition  process  in¬ 
variably  lead  to  chain-branching  processes  by  for¬ 
mation  either  of  OH  and  a  carbonylhydroperoxide¬ 
like  species  (reaction  7)  or  of  HO2  and  an  olefinic 
hydroperoxide  species.  The  stable  species  formed  in 
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Fig.  5.  Key  low-  and  high-temperature  reaction  path¬ 
ways  of  methyl  butanoate  for  one  hydrogen  abstraction  site. 
Carbon  atoms  are  assumed  to  be  saturated  with  H  atoms, 
except  where  radical  site  is  noted. 


this  way  readily  decomposes  to  form  two  radicals 
(reaction  8). 

In  each  of  these  competitions,  the  rate  of  the  uni- 
molecular  decomposition  process  (producing  no  net 
radicals)  increases  more  rapidly  with  increasing  tem¬ 
perature  than  that  of  the  O2  addition  process  (ulti¬ 
mately  producing  as  many  as  two  net  radicals).  Thus, 
although  all  rates  increase  as  temperature  increases, 
higher  temperatures  favor  the  chain-propagating 
processes  over  the  chain  branching  ones,  leading  to 
a  net  decrease  in  reactivity  over  some  temperature 
range.  This  change  can  be  illustrated  for  the  condi¬ 
tions  of  Figs.  1  and  3  by  considering  the  fate  of  the 
radical  formed  by  abstraction  of  one  of  the  H  atoms 
beta  to  the  carbonyl  group  in  methyl  butanoate  (Fig. 
5).  At  665  K,  at  the  point  where  half  the  fuel  is  con¬ 
sumed,  over  99%  of  this  species  reacts  with  O2  (re¬ 
action  3).  At  765  K,  unimolecular  decomposition  to 
methyl  and  an  olefinic  methyl  ester  (reaction  2)  be¬ 
comes  the  dominant  route,  accounting  for  75%  of 
the  reaction  flux,  and  at  865  K  the  unimolecular  de¬ 
composition  pathway  accounts  for  more  than  97%  of 
the  reaction  flux.  The  fate  of  the  isomerization  prod¬ 
ucts  of  the  initial  O2  addition  reaction  changes  with 
temperature  as  well.  At  665  K,  over  50%  of  the  flux 
through  these  species  leads  to  chain  branching  (re¬ 
action  6),  while  at  765  K  only  25%  leads  to  chain 


TABLE  2 

Arrhenius  Parameters  (k  =  AT"  exp{  —  E/^/{RT)))  and  Sensitivity  Coefficients  of  Selected  Reactions 


Reaction 

Sensitivity 
Coefficient 
at  665  K 

Sensitivity 
Coefficient 
at  765  K 

A 

(cm^  mol,  s  units) 

n 

(for  Tin  K) 

Ea 

(cal/mol) 

mb  +  o2  <=>  ho2  +  mb2j^ 

0.083 

0.055 

4.OOE  4-  13 

0 

41,300 

mb2oo  mb2ooh4j 

-  0.068 

-0.011 

6.10E  4-  10 

0 

25,076 

mb  +  lio2  <=>  h2o2  +  mb2j 

-  0.054 

-  0.059 

4.32E  4-  12 

0 

14,400 

mb  +  oh  ^  h2o  +  mb3j 

0.048 

0.020 

4.68E  4-  07 

1.61 

-35 

mb  +  oh  h2o  +  mb4j 

-  0.044 

-  0.007 

5.25E  4-  09 

0.97 

1,590 

mb2oo  mb2ooh3j 

0.041 

0.016 

1.09E  4-  12 

0 

29,206 

mb  +  ch3o2  ^  ch3o2h  +  mb2j 

-  0.028 

-  0.055 

4.00E  4-  12 

0 

14,000 

mb  +  mb2oo  <=>  mb2ooh  +  mb2j 

-  0.026 

abs  value  <0.007 

2.16E  4-  12 

0 

14,400 

ch3o  +  o2  ^  ch2o  +  ho2^ 

0.023 

0.024 

5..50E  4-  10 

0 

2,424 

co2  +  ch3  ^  ch3oco^ 

-  0.020 

abs  value  <0.007 

1..50E  4-  11 

0 

36,730 

mp2d  +  ch3  ^  mb2j^ 

abs  value  <0.007 

-  0.052 

l.OOE  4-  11 

0 

7,600 

ch3o2  +  ch3  <=>  2ch3o^ 

abs  value  <0.007 

0.016 

7.00E  4-  12 

0 

- 1,000 

mb  +  oh  ^  h2o  +  mb2j 

-0.015 

-  0.015 

1.15E  4-  11 

0.51 

63 

mb  +  ch3o  <=>  ch3oh  +  mb2j 

-0.013 

-  0.013 

3.80E  4-  10 

0 

2,800 

Note:  Sensitivity  coefficients  represent  the  sensitivity  of  the  methyl  butanoate  concentration  to  change  in  individual 
rate  constant  using  standard  CHEMKIN  sensitivity  analysis  [14],  at  the  point  where  10%  of  the  fuel  is  consumed. 
Conditions  the  same  as  Fig.  1.  Initial  temperatures  are  665  and  765  K.  “mb”  denotes  methyl  butanoate  structure,  “mp” 
denotes  methyl  propanoate  structure,  “j”  denotes  radical  site.  Numbered  label  (e.g.,  mb2j)  refers  to  carbon  position  where 
the  C  =  0  carbon  is  carbon  1.  “d”  denotes  double  bond.  For  hydroperoxyalkyl  radical  analogs  (e.g.,  mb2ooh4j),  the  first 
number  refers  to  “ooh”  location  and  the  second  number  to  radical  site  location. 

“reaction  proceeds  in  reverse  of  direction  written  here. 

'’reverse  rate  constants  specified  as  well. 
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branching.  Methyl  formate  follows  a  similar  pattern, 
with  a  shift  from  reaction  with  O2  addition  and  chain 
branching  to  chain-propagating  unimolecular  de¬ 
composition  reactions  when  temperatures  increase 
through  the  NTC  region. 

Table  2  shows  the  reactions  to  which  methyl  bu- 
tanoate  is  most  sensitive  at  temperatures  in  and 
above  the  NTC  region.  Both  sensitivity  coefficients 
and  Arrhenius  parameters  are  listed.  Reactions  in¬ 
volving  the  carbon  beta  to  the  carbonyl  group  are 
prominent  at  both  temperatures  because  of  the  rela¬ 
tively  weak  C-H  bond  there.  A  surprising  result  is 
the  sensitivity  to  the  reaction  of  HO2  with  a  key 
methyl  butanoate  radical,  a  termination  reaction 
which  slows  down  the  overall  rate  of  oxidation.  Isom¬ 
erization  reactions  such  as  reaction  4  (Fig.  5)  are 
important  in  the  NTC  region  (665  K).  Above  the 
NTC  region  (765  K),  reactions  of  the  fuel  with  RO2 
radicals  control  the  chemistry  and  accelerate  the 
overall  rate  of  oxidation. 


Summary  and  Conclusions 

Detailed  chemical  kinetic  models  for  combustion 
of  methyl  butanoate  and  methyl  formate  have  been 
developed  following  rules  developed  for  hydrocar¬ 
bon  combustion.  The  mechanisms  have  been  tested 
against  the  limited  combustion  data  available  in  the 
literature,  which  was  obtained  at  low  temperature, 
subatmospheric  conditions  in  static  vessels,  using 
pressure  measurements  as  the  main  diagnostic. 
Some  qualitative  agreement  was  obtained,  but  the 
experimental  data  consistently  indicated  lower  over¬ 
all  reactivities  than  the  model,  differing  by  factors  of 
10  to  50.  This  discrepancy,  which  occurs  for  species 
with  well-established  kinetic  mechanisms  as  well  as 
for  methyl  esters,  is  tentatively  ascribed  to  the  pres¬ 
ence  of  wall  reactions  in  the  experiments.  The  model 
predicts  a  region  of  weak  or  negative  dependence  of 
overall  reaction  rate  on  temperature  for  each  methyl 
ester.  Examination  of  the  reaction  fluxes  provides  an 
explanation  of  this  behavior,  involving  a  tempera¬ 
ture-dependent  competition  between  chain-propa¬ 
gating  unimolecular  decomposition  processes  and 
chain-branching  processes,  similar  to  that  accepted 
for  hydrocarbons.  There  is  an  urgent  need  to  obtain 
more  complete  and  well-characterized  experimental 
data  for  thorough  testing  of  the  model. 
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COMMENTS 


William  H.  Green,  MIT,  USA.  Liquid-phase  low-tem¬ 
perature  oxidation  experiments  indicate  that  the  C-H  bond 
adjacent  to  a  carbonyl  group  in  an  ester  is  stronger,  not 
weaker,  than  an  alkyl  C-H.  It  may  be  worthwhile  to  reex¬ 
amine  the  thermochemical  calculations  on  this  bond 
strength. 

Author’s  Reply.  We  have  two  computed  values  of  the 
bond  strength  available.  One  is  96.2  kcal/mol  based  indi¬ 
rectly  on  BAC  MP4  calculations  of  Melius  (Ref  [20]  in  this 
paper)  as  explained  in  the  paper  (Table  1);  the  other  is  93.1 
kcal/mol  calculated  by  Rozzelli  and  Chen  (Ref  [18]  in  this 
paper)  using  the  CBS-Q  method  with  isodesmic  reactions. 
Therefore,  we  think  it  is  unlikely  that  the  bond  strength  is 


stronger  than  a  secondary  alkyl  C-H  (98.5  kcal/mol)  (Ref 
[17]  in  this  paper).  It  would  be  useful  to  have  more  details 
on  the  nature  of  the  experimental  evidence  that  you  men¬ 
tion  for  a  stronger  C-H  bond.  Presumably,  the  experiments 
rely  on  measurements  of  stable  products,  and  their  results 
may  be  attributable  to  kinetic  or  other  factors  besides  dif¬ 
ferences  in  bond  strength.  However,  a  reexamination  of  the 
calculations  would  be  worthwhile  to  see  if  some  detail  was 
overlooked.  Our  understanding  of  these  liquid-phase  ex¬ 
periments  is  that  the  conclusion  of  “higher”  or  “lower” 
bond  strength  is  inferred  from  product  analysis  of  the  ox¬ 
idation  of  the  ester.  It  may  also  be  worthwhile  to  reexamine 
the  interpretation  of  these  experiments  that  give  an  indirect 
indication  of  relative  bond  strength. 


